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Abstract 



The ALICE Collaboration has measured inclusive J/\j/ production in pp collisions at a center of mass 
energy y/s = 2.76 TeV at the LHC The results presented in this Letter refer to the rapidity ranges 
\y\ < 0.9 and 2.5 < y < 4 and have been obtained by measuring the electron and muon pair decay 
channels, respectively. The integrated luminosities for the two channels are L? t = 1 . 1 nb~ 1 and L^ m — 

19.9 nb _1 , and the corresponding signal statistics are Ny*~ = 59 ± 14 and Nj, J 1 = 1364±53. We 

present dOj/^/dy for the two rapidity regions understudy and, for the forward -grange, d 2 (7j j y / dydp t 
in the transverse momentum domain < p x < 8 GeV/c. The results are compared with previously 
published results at *Js = 7 TeV and with theoretical calculations. 
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1 Introduction 

Almost forty years after the discovery of charmonium, its production in hadronic collisions still remains 
not completely understood, and charmonium production data represent a complex and severe test for 
QCD-inspired models HI. 

Recently, first results from the Large Hadron Collider (LHC) on J/i/a production in pp collisions at y/s = 
7 TeV became available |[2j-[6), significantly extending the energy reach beyond that of the Tevatron and 
RHIC hadron colliders [7H9]]. A reasonable description of the transverse momentum spectra has been 
obtained by theoretical models ifTOHTBl . and first results on J/ y polarization, a crucial testing ground for 
theory lfl4l - [T6Tl . are also available ifTTl at LHC energy. 

At the beginning of 2011, the LHC delivered pp collisions at y/s = 2.76 TeV. The main goal of this short 
run was to provide a reference for the Pb-Pb data which were taken at the same yfs per nucleon-nucleon 
collision. On the other hand, these data offer the possibility of studying J/y production at an intermediate 
energy between Tevatron and the present LHC top energy, and represent therefore an interesting test for 
models. 

In this Letter, we present results on inclusive J/y production at y/s = 2.76 TeV as obtained by the 
ALICE experiment iTTBl . J/y particles were measured, down to zero transverse momentum, via their 
decay into e + e~ at mid-rapidity (|y| < 0.9) and into jU + jU~ at forward rapidity (2.5 < y < 4). Results 
from ALICE on J/y production at y/s = 7 TeV were recently published OUT]. Since the experimental 
apparatus and the data analysis procedure are basically the same for the two data samples, they will be 
concisely described, referring where necessary to our previous publications. Results will then be shown 
for dCj/^/dy at central and at forward rapidity. For the region 2.5 < y < 4 the differential cross section 
d 2 Oy v /dydp t will also be given, for the transverse momentum range < p t < 8 GeV/c. A comparison 
with the previous results at y/s = 7 TeV will be carried out and next-to-leading order Non-Relativistic 
QCD (NLO NRQCD) theoretical calculations will be compared to the experimental data. 

2 Experimental apparatus and data analysis 

The main elements of the ALICE experiment at the CERN LHC are a central rapidity barrel (covering 
the pseudorapidity range |tj| < 0.9) for the detection of hadrons, electrons and photons and for the 
measurement of jets, and a forward muon spectrometer (—4 < r\ < —2.5). The experimental set-up is 
described in detail in iTTBl . For the analysis described in this Letter, the relevant detector systems for 
tracking and electron identification in the central barrel are the Inner Tracking System (ITS) |fl9l , based 
on six layers of silicon detectors, and the Time Projection Chamber (TPC) BUI . The ITS covers the 
|t] | < 0.9 range and, together with two small forward scintillator detectors (VZERO, covering 2.8 < r\ < 
5.1 and —3.7 < r\ < —1.7), is used to define the Minimum-Bias (MB) interaction trigger. In particular, 
the MB condition requires a logical OR between at least one fired read-out chip in one of the two inner 
layers of the ITS (Silicon Pixel Detector), and a signal in at least one of the VZERO detectors. Muons are 
tracked and identified in the muon spectrometer [5 1, which consists of a front absorber to remove hadrons, 
a 3 Tm dipole magnet and a tracking system. It also includes a triggering system with a programmable 
p t threshold. With this trigger, the collected data sample was enriched with events where, in addition 
to the MB condition, at least one muon was detected in the spectrometer acceptance. The threshold for 
the muon trigger was set to its minimum value, p t = 0.5 GeV/c. With this choice the acceptance for 
J/y — > £i + jU~ detection extends down to p t = 0. Further details on the detectors relevant for this analysis 
and on the trigger definitions can be found in Ref. 0. 

The dielectron analysis is based on a sample of 65.4- 10 6 MB triggers, corresponding to an integrated 
luminosity L? nt = 1.1 nb _1 . Out of the total sample, 47.4-10 6 events have a reconstructed vertex which 
lies within ±10 cm, along the beam direction, from the nominal interaction point and are retained for 
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the following analysis steps. The analysis strategy is briefly described below. It is the same as applied in 
case of the analysis at y/s = 7 TeV, small differences are explained in the text. For details we refer to 0. 

Reconstructed tracks are required to have a hit in one of the two innermost or in the fifth ITS layer (layers 
three and four were excluded from the reconstruction). This choice makes the track cuts somewhat less 
stringent as compared to the analysis of the y/s = 7 TeV data where a hit was required in one of the 
two innermost layers. As a result, the signal increases by ~12%, whereas the significance for the two 
cuts is comparable within the uncertainties. The choice to use the looser cut was motivated by the 
fact that it provides a central cross section value of the systematic variations using different cuts. The 
number of TPC clusters for each track must be larger than 70 (out of a maximum of 159), with the % 2 
per space point of the momentum fit lower than 4. The kinematic cuts p t > 1 GeV/c and |t]| < 0.9 
are applied to each track. The electron identification is based on the correlation between the specific 
energy loss dE/dx and the momentum measured in the TPC, requiring a ±3a inclusion cut around the 
electron line corresponding to the Bethe-Bloch expectation and an exclusion cut of ±3. 5a (±3(7) for 
pions (protons). Finally, a rapidity cut \y\ < 0.9 is applied to to J/i/a candidates to remove pairs at the 
edge of the acceptance. 

The signal extraction is based on the like-sign (LS) subtracted invariant mass spectrum of e + e~ pairs. The 
LS spectrum is obtained as the sum of positive-positive and negative-negative spectra. The scale factor 
on the LS background, applied in Q to account for various non-combinatorial effects, was found to be 
negligible in this analysis. Figure Q] (top panel) shows the opposite sign (OS) dielectron mass spectrum 
together with the LS spectrum. After subtraction, the number of J/y is estimated by bin counting in the 
invariant mass range 2.92 < m e + e - < 3.20GeV/c 2 , resulting in 59 ± 14 (stat.) counts with a significance 
of 5.4 ± 0.6. The signal fraction in the mass range defined above is estimated from a Monte Carlo 
(MC) simulation, and included in the acceptance. In Fig. Q] (bottom panel) the LS-subtracted spectrum 
is overlayed with the MC signal shape, normalized to the data points in the invariant mass range 2.5 < 
m e + e - < 3.5GeV/c . In addition to the LS method, the background estimated using a track rotation 
(TrkRot) technique^ is also shown in Fig. Q] The differences between the number of J/yf obtained with 
the TrkRot and LS methods is used in the estimate of the systematic uncertainty on the signal extraction. 

The dimuon analysis is based on 8.8T0 6 muon-triggered events, corresponding to an integrated lumi- 
nosity L^ t = 19.9 nb _1 . Out of this sample, 1.0T0 5 events contain a reconstructed OS muon pair. It 
is required that each event contains at least one reconstructed vertex. Events are retained for the anal- 
ysis if both candidate muon tracks exit the front hadron absorber (z = —503 cm) at a radial coordinate 
17.6 < 7? a bs < 89.5 cm, a cut roughly corresponding to the angular acceptance of the muon spectrometer. 
It is also required that at least one of the two muons satisfies the muon trigger condition. Finally, the cut 
2.5 < y < 4 is applied to the pairs in order to reject dimuons at the edge of the spectrometer acceptance. 

The signal is extracted by a fit to the invariant mass spectrum over the range 2 < < 5 GeV/c . The 
signal is parameterized with a Crystal Ball (CB) function [21] with a background described by the sum 
of two exponentials. The position (m 3 /y) of the peak of the CB function, as well as its width (wj/^), are 
kept as free parameters in the fit. The obtained values are my v = 3.129± 0.004 GeV/c 2 (a value larger 
by ~1% than the world average J22j]) and w J/v = 0.083 ±0.004 GeV/c 2 . The J/i/a width is only slightly 
larger (by ~ 0.006 GeV/c 2 ) than that obtained in the MC, which includes the effect of the misalignment 
of the muon tracking system. The tails of the CB function are fixed to their MC value, since with 
the available statistics and signal to background ratio they cannot be reliably extracted from the fitting 
procedure. Finally, the contribution of the y{2S) signal is included in the fit, although its influence on the 
number of detected J/i/a is negligible. In Fig.[2]the dimuon invariant mass spectrum is presented, together 
with the result of the fit {% 2 jndf = 1.3). By integrating the CB function, one gets a total number of J/i/a 



In the TrkRot method one track of the OS pair is rotated around the z-axis. The procedure is repeated several times 
randomly varying the rotation angle. In this way,one removes the correlation between the two electrons of the pair. The TrkRot 
invariant mass spectrum is scaled to match the integral of the OS spectrum in the region 3.2 < m e + e - < 5 GeV/c 2 . 
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Fig. 1: Top panel: invariant mass distributions for opposite-sign (OS) and like-sign (LS) electron pairs (\y\ < 0.9, 
all px). The background estimate from the TrkRot method (see text for details) is also shown. Bottom panel: the 
difference of the OS and LS distributions with the normalized MC signal shape superimposed. 



N?,» = 1364 ±53 (stat). 

The J/y statistics in the dimuon channel permit a differential study of the production cross sections 
using six y or seven p t intervals. The fitting technique is the same as for the integrated invariant mass 
spectrum, except for the value of the CB width which was fixed for each bin i to the value Wj,^ = 

w J/v ' ( w j/> : / w j/v/)' i- e -' by scaling the measured width for the integrated spectrum with the MC ratio 
between the widths for the bin i and for the integrated spectrum. The sum of the signal events for both p t 
and y bins agrees well (within 0.3% and 1.2% respectively) with the result of the fit to the integrated mass 
spectrum. In Fig. [3] the invariant mass spectra corresponding to the various p t bins are shown, together 
with the results of the fits. The J/y signal is well visible also in the spectra with lower statistics and the 
quality of the fits is similar to the one obtained for the integrated mass spectrum. 

For both the dielectron and dimuon analyses the number of signal events is corrected by the product of 
acceptance times efficiency (A x e). The Axe values are obtained using MC simulations which include 
a description of the status of the detector as a function of time. Details on the procedure are given in 
Ref. Q. For this analysis, the MC input distributions in transverse momentum and rapidity are obtained 
by interpolating between the LHC results for yfs = 7 TeV and lower energy collider measurements ll23l . 
It was verified a posteriori that the interpolated input spectra are in good agreement with those obtained 
from this analysis. The results are A x e = 0.120 for the dielectron analysis and A x e = 0.346 for the 
dimuon analysis. These values refer to J/y production for p t > in the analyzed rapidity ranges, |y| < 0.9 
and 2.5 < y < 4, respectively. 

The inclusive J/y production cross section for the leptonic channel £ + £~ is calculated as: 
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Fig. 2: Invariant mass distribution for opposite-sign muon pairs (2.5 < y < 4, all p t X in the mass region 2 < mnu < 5 
GeV/c 2 , with the result of the fit (see text for details). The fitted J/y/ and y(2S) contributions, as well as the 
background, are also shown. 
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where A 7 ?^ ^ = A 7 -^ /(A x £) f+( is the number of signal events corrected for acceptance times effi- 
ciency, BR(J/y^ l + £~) = (5. 94 ±0.06)% E3 is the leptonic branching ratio for the J/y decay, N M b is 
the number of MB -triggered events and (7mb = 55.4 ± l.O(total) mb is the absolute cross section for the 
occurrence of the MB condition [24], derived from the result of a van der Meer scan (see [5 ] for details). 
The R (+e factor is 1 for the e + e~ analysis, whereas for the dimuon channel R^ + ^ = 0.0326 ± 0.0002 
represents the inverse of the enhancement factor of the muon trigger with respect to the MB trigger |0. 
An equivalent formula is used for the differential cross sections in y and p t . 

The sources of systematic uncertainties are exactly the same as for the corresponding y/s = 7 TeV anal- 
ysis and have been estimated in a similar way (see Q for details). In Table [T] we quote their values 
for the integrated cross sections in the dielectron and in the dimuon channel. The uncertainty on sig- 
nal extraction for the electron analysis (14%) is larger than the 8.5% quoted at y/s = 7 TeV Q. This 
increase mainly comes from the difference in A 7 ?^ 6 e obtained by requiring various conditions in the 
ITS: a hit in the first layer, in any of the first two layers (as was done for the = 7 TeV analysis), or 
the less stringent condition adopted from the present analysis, described earlier in this Section. For the 
muon analysis, the uncertainty on signal extraction (4%) is now smaller with respect to the 7.5% quoted 
at y/s = 7 TeV 0. The present value was calculated as the average absolute deviation on the number 
of signal events obtained with alternative parameterizations of the signal and background shapes. At 
y/s = 1 TeV the more conservative, but also more prone to statistical effects, approach of using the larger 
deviation obtained in the various fits was adopted. Finally, the decrease of the systematic uncertainty on 
the trigger efficiency for the muon analysis (from 4% at y/s = 7 TeV Q to 2% at y/s = 2.76 TeV) is due 
to a different approach, the present one being based on the study of the variation of the J/y/ triggering 
efficiency when the efficiency of the trigger detectors is changed by an amount slightly larger than the 
uncertainty on this last quantity. 

The total systematic uncertainties, excluding those related to the unknown degree of polarization of the 
J/y/", are 18.0% and 8.1% for the dielectron and the dimuon channel, respectively. For the differential 
cross sections measured in the dimuon channel, the same sources of systematic uncertainties quoted in 
Table [Q apply to each y and p t bin. For the uncertainties relative to the choice of the MC inputs, their 
values may in principle vary with either rapidity or transverse momentum. However, no clear trend as 
a function of these two variables is observed. So, the relative systematic uncertainty calculated for the 
integrated cross section is assigned to each point and considered as uncorrelated between the bins. The 
uncertainty on signal extraction is also considered as bin-to-bin uncorrelated. The limited signal statistics 
for most of the bins prevents a direct study of the systematic uncertainty, therefore the relative systematic 
uncertainty assigned to the integrated cross section was attributed to each point. 



3 Results 

The analysis described in the previous Section gives the following results for the integrated inclusive J/\j/ 
cross sections in the two rapidity ranges investigated at y/s = 2.76 TeV: 

Oj/y(M < °- 9 ) = 7 - 75 ± l-78(stat.) ± 1.39(syst.) + 1.16(A HE = 1) - 1-63(A HE = -1) fib and 
<Jj /v (2.5 <y<4) = 3.34 ± 0.13(stat.) ± 0.27(syst.) + 0.53(Acs = 1) - 1-07 (Acs = -1) jub. 

The polarization-related systematic uncertainties were estimated in the helicity (HE) and Collins-Soper 
(CS) reference frames [25 1 . The uncertainties are quoted in the frames where they are larger. Existing 
polarization results for yfs = 7 TeV at forward rapidity [17 1, tend to exclude a significant degree of 
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Table 1: Systematic uncertainties (in percent) contributing to the measurement of the integrated J/y/ cross section. 
The uncertainties related to the J/y/ polarization were calculated for both Collins-Soper and helicity reference 
frames. 



Channel 


e + e 




Signal extraction 


14 


4 


Acceptance input 


1.5 


4 


Trigger efficiency 




2 


Reconstruction efficiency 


11 


4 


R factor 




3 


Luminosity 


1.9 


1.9 


B.R. 


1 


Polarization 


A = -1 


X = \ 


X = -\ 


X = 1 


CS 


+19 


-13 


+32 


-16 


HE 


+21 


-15 


+24 


-12 



polarization for the J/y/. However, in absence of clear predictions for the y^-dependence of the effect, 
we prefer to quote systematic uncertainties relative to fully longitudinal (X = — 1) or transverse (X = 1) 
degree of polarization. With respect to the y/s = 7 TeV measurement, the y/s = 2.76 TeV cross sections 
are smaller by a factor 1.59±0.50 (1.89±0.31) for the \y\ < 0.9 (2.5 < y < 4) rapidity ranges. The 
quoted uncertainty on the ratios is obtained by propagating the quadratic sum of statistical and systematic 
uncertainties (excluding the polarization-related contribution) of the two cross section values. 

Figure |4] presents the differential cross section d 2 Gj/ Xj/ /dp t dy, averaged over the interval 2.5 < y < 4, 
for the transverse momentum range < p t < 8 GeV/c. The results are compared with those previously 
published by ALICE for y/s = 7 TeV, as well as, for the range 3 < p t < 8 GeV/c, with the predictions 
of a NRQCD calculation 1251 . which includes both colour singlet and colour octet terms at NLO. The 
model satisfactorily describes both sets of experimental data. 

Using the results shown in Fig.|4l the mean transverse momentum for inclusive J/ty/ production at forward 
rapidity is computed by fitting d 2 Oj/ ¥ /dp t dy with the function 



" :CT -C , » .... (2) 



dptdy 



l + ( & 

1 Po 



with C, po and n as free parameters, as done in ||9]. The result, relative to the range < p t < 8 GeV/c, 
is (p t ) = 2.28 ±0.07(stat) ±0.04(syst) GeV/c. A similar analysis carried out on the y/s = 7 TeV data 
published in Q gives (p t ) = 2.44 ± 0.09(stat) ± 0.06(syst) GeV/c for 2.5 < y < 4 and (p t ) = 2.72 ± 
0.21(stat) ±0.28(syst) GeV/c for |y| < 0.9 (for that data sample d 2 Oyy/dp t dy was also calculated for 
the dielectron analysis, in the range < p t < 7 GeV/c). The quoted systematic uncertainties are related 
to the uncorrelated systematic uncertainties for d 2 a /dp t dy. 

Figure [5] presents the y^-dependence of the inclusive J/y (pt), for various fixed-target and collider 
experiments |[3ll5Tr71l9ll27l. The results show a roughly linear increase of (p t ) with ]n(y/s), with slightly 
larger (p t ) values at central rapidity. The numerical values for both (p t ) and (p 2 ) are quoted in Table[2] 

In Fig.[6]we present the results for do^y/dy at = 2.76 TeV, compared with the previously published 
^/s = 7 TeV results. The numerical values corresponding to the results presented in Fig.|4]and Fig.[6]are 
shown in Table [3j together with the number of signal events and with the values for A x e. Most sources 
of systematic uncertainty are common or strongly bin-to-bin correlated, except, as outlined before, the 
ones related to the signal extraction and to the MC inputs that are therefore quoted separately in Table [3] 
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Fig. 4: Double differential J/y/ production cross section at y/s = 2.76 TeV compared to previous ALICE results at 
^/s = 7 TeV [5 1. The vertical error bars represent the statistical errors while the boxes correspond to the systematic 
uncertainties. The systematic uncertainties on luminosity are not included. The results are compared with a NLO 
NRQCD calculation 1 26 1 performed in the region p t > 3 GeV/c. 
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Fig. 5: The -^-dependence of (p t ) for inclusive J/\j/ production, for various fixed-target and collider experiments. 
For the ALICE points the error bars represent the quadratic sum of statistical and systematic uncertainties. The 
points for y/s = 7 TeV have been slightly shifted to improve visibility. 
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Table 2: The (p t ) and {pi) values for inclusive i/y/ production measured by ALICE. Statistical and systematic 
uncertainties are quoted separately. 





(pt) (GeV/c) 


04) (GeV/c) 2 


y/s = 2.76 TeV, 2.5 < y < 4 


2.28 ± 0.07 ± 0.04 


7.06 ± 0.40 ± 0.22 


y/s = 7 TeV, |y| < 0.9 


2.72 ± 0.21 ± 0.28 


10.02 ± 1.40 ± 1.80 


yfi = 7 TeV, 2.5 < y < 4 


2.44 ± 0.09 ± 0.06 


8.32 ± 0.50 ± 0.35 



T3 
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Fig. 6: Differential ]ly production cross section at y/s = 2.76 TeV compared to previous ALICE results at y/s = 
7 TeV [5|. The vertical error bars represent the statistical errors while the boxes correspond to the systematic 
uncertainties. The systematic uncertainties on luminosity are not included. 



Table 3: Summary of the results concerning the J/y/ differential cross sections for pp at s/s = 2.76 TeV. 



Pt 




Axe 


d 2 a ]/v /dp t dy 




Systematic uncertainties 




(GeV/c) 






Ob/(GeV/c)) 


Correl. 


Non-correl. 


Polariz., CS 


Polariz., HE 










(Aib/(GeV/c)) 


(/ib/(GeV/c)) 


(/ita/(GeV/c)) 


(Mb/(GeV/c)) 


2.5<y<4 


[0;i] 


222±19 


0.330 


0.380±0.033 


0.022 


0.021 


±0.074,-0.141 


±0.069,-0.133 


[i;2] 


407±24 


0.326 


0.705 ±0.042 


0.041 


0.040 


±0.122,-0.271 


±0.098,-0.211 


[2; 3] 


343±22 


0.332 


0.583±0.038 


0.034 


0.033 


±0.100,-0.203 


±0.069,-0.127 


P;4] 


201±17 


0.354 


0.321 ±0.027 


0.019 


0.018 


±0.050,-0.089 


±0.029, -0.047 


K;5] 


95±12 


0.397 


0.135±0.017 


0.008 


0.008 


±0.014,-0.027 


±0.009,-0.018 


[5;6] 


58±9 


0.449 


0.073±0.011 


0.004 


0.004 


±0.005,-0.011 


±0.005,-0.009 


[6; 8] 


34±7 


0.507 


0.019±0.004 


0.001 


0.001 


±0.001,-0.001 


±0.001,-0.002 


y 






da ]/¥ /dy (jib) 


(Mb) 


(lib) 


(Mb) 


(Mb) 


[-0.9;0.9] 


59±14 


0.120 


4.31±0.99 


0.08 


0.77 


±0.57,-0.81 


±0.65,-0.90 


[2.5;2.75] 


121±14 


0.134 


3.05±0.35 


0.18 


0.17 


±0.67,-1.41 


±0.52,-1.04 


[2.75; 3] 


252±20 


0.361 


2.37±0.19 


0.14 


0.13 


±0.42,-0.84 


±0.39,-0.78 


[3; 3.25] 


325±22 


0.488 


2.26±0.15 


0.13 


0.13 


±0.29,-0.65 


±0.31,-0.61 


[3.25;3.5] 


298±21 


0.502 


2.01±0.14 


0.12 


0.11 


±0.27,-0.54 


±0.21,-0.38 


[3.5; 3.75] 


245±19 


0.416 


2.00±0.16 


0.12 


0.11 


±0.33,-0.67 


±0.15,-0.30 


[3.75; 4] 


106±12 


0.214 


1.68±0.19 


0.10 


0.09 


±0.36,-0.69 


±0.16,-0.26 
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The kinematic coverage of the ALICE experiment is unique among the LHC experiments due to the very 
good acceptance down to p t = at central rapidity. This feature allows a comparison of the /^-integrated 
midrapidity cross sections with those from lower energy collider experiments. The result is displayed in 
Fig. [TJ where the dCj/^/dy values from ALICE for the two energies are shown together with results from 
RHIC [9] and Tevatron [7] experiments, as a function of y/s. 
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Fig. 7: -^-dependence of the inclusive }/y/ production cross section da/dy, at central rapidity for various collider 
experiments. 



4 Conclusions 

The ALICE experiment has measured the inclusive J/i/a production cross section for proton-proton colli- 
sions at y/s = 2.76 TeV, in the rapidity ranges \y\ < 0.9 and 2.5 < y < 4, down to p t = 0. The measured 
values are o ]/v (\y\ < 0.9) = 7.75 ± 1.78(stat) ± 1.39(syst) + 1.16(A HE = 1) - L63(A HE = -1) lib and 
a J/v/ (2.5 <y < 4) = 3.34 ±0.13(stat) ±0.27(syst) + 0.53(Acs = 1) - 1.07 (Acs = -1) jub. Differen- 
tial cross sections in y and p t have also been measured for the forward rapidity region. These results 
provide an important intermediate point between top Tevatron energy and the current maximum LHC 
energy. They also represent a crucial reference for the measurement of nuclear effects on J/y/ production 
in Pb-Pb interactions carried out at the same centre-of-mass energy per nucleon-nucleon collision [28 1 . 
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